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SOIL-PILE INTERACTION IN HORIZONTAL VIBRATION 

M. NOVAK* AND T. NOGAMlf 
Faculty of Engineering Science, The University of Western Ontario, London, Canada 

SUMMARY 

Interaction between soil and an elastic pile vibrating horizontally is theoretically examined. The soil is modelled as a 
linear, viscoelastic layer overlying rigid bedrock. The pile is assumed to be vertical and point bearing. This study 
utilizes the definition of soil resistance presented in a preceding paper. 1 A direct solution is developed which yields 
closed form formulas for pile displacement, stiffness and damping. A parametric study clarifies the role of the 
parameters involved, illustrates the interaction between the soil and the pile and shows the stiffness and damping 
properties of the soil-pile system for typical values of the governing parameters. 

INTRODUCTION 

A pile exposed to dynamic horizontal loads, or eccentric vertical loads, such as those due to machines, wind 
or earthquakes, interacts with the surrounding soil. The result of this interaction is the modification of pile 
stiffness and generation of damping through energy radiation and dissipation. Energy radiation results from 
the propagation of progressive waves. Energy dissipation is caused by material damping of the soil and, to a 
lesser degree, of the pile and by friction between the pile and the soil. 

In this paper, the modification of pile stiffness and the generation of damping due to energy radiation and 
dissipation in the soil are analysed and the response of the pile to horizontal loads and moments in the vertical 
plane is theoretically investigated. This study is based on the definition of the resistance of the soil layer to 
the pile motion presented in a preceding paper. 1 The main assumptions of the solution are: 

The soil layer is linear, viscoelastic, free at the surface and fixed at the bottom. The material damping is of 
the hysteretic frequency independent type. The pile is vertical, uniform, linearly elastic and of circular cross- 
section. It is clamped or pinned at the tip and perfectly attached to the soil 

The same assumptions were adopted in Reference 2 when analysing the vertical response. Similar assump- 
tions were made by Tajimi 3 who assumed frequency dependent viscous damping (the Kelvin-Voigt model) 
and solved the horizontal response of a pile assuming a fixed head (no rotation) and a pinned tip. Novak 4 
assumed a simpler model in which only horizontally propagating waves were considered. This assumption is 
identical to that of plane strain conditions; it yields simple formulas and good accuracy for slender piles and 
for frequencies higher than the first natural frequency of the soil layer. Penzien 5 employed a discretized 
model which made it possible to incorporate non-linearity as this appears to be desirable for large amplitudes. 

The advantage of the approach presented is that it offers basic insight into the mechanism of the soil-pile 
interaction. The results can be used for guidance in design or when choosing parameters for discretized 
models. 

VIBRATION OF THE SOIL-PILE SYSTEM 

When the pile vibrates harmonically due to external excitation it meets the resistance of the surrounding soil 
(Figure 1). For the assumptions adopted in this paper the soil resistance was defined in Reference 1 as a 
function of a given motion of the pile. This motion is not known beforehand but it can be found from the 
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condition that the motion of the soil at the pile surface, i.e. at stations (r 0 ,z), must be equal to the motion 
of the pile. Two approaches were used to solve this problem: modal analysis and a direct solution. The modal 
analysis can be readily employed using the mode shapes of a free-standing pile and considering the soil 
resistance as an external load. This approach was used by Tajimi 8 and initially also by the writers. The direct 
method used in this paper is much more efficient than modal analysis in that it always leads to small matrices, 
the size of which is only 4x4. Consequently, the computing costs are low. 
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Figure 1. External forces and soil resistance to motion of pile 



Assume that the pile is subject to harmonic excitation by end forces applied at the head of the pile (Figure 
1). Then the pile motion is induced through the boundary conditions. If the horizontal motion of the pile is 

u(z,t) = u(z)eM (!) 
where u = u(z) is the complex amplitude, the governing equation of the pile motion is 

E p I^(u^)+m~(u^) = -/>(z)e^ (2) 

in which E p I- bending stiffness of the pile, m = mass of the pile per unit length and p(z) « amplitude of the 
soil resistance to the motion of the pile. Under the assumptions adopted, the soil resistance is given by equa- 
tion (26) in Reference 1 and is 

00 

P(z)= 2 oc hn U n smh n z (3) 

where a hn = horizontal resistance factor depending on pile radius r 0 , shear modulus fx and a number of 
dimensionless parameters such as the dimensionless frequency a 0 = Ha>lv 89 pile slenderness Hjr 09 material 
damping D and Poisson's ratio v; U n = modal amplitude independent of z; sin h n z = mode shapes of the 

soil layer and h n = (7r/2H)(2n-l) where H = depth of the layer and = 1, 2, 3, 

Substitution of equation (3) into equation (2) and elimination of time t yields the following equation for the 
pile amplitude, 

d 4 w 0 00 
E v I 'J^~ ma) " = ~ S oc hn U n sinh n z (4) 

71=1 

The solution to equation (4) is a sum of the complete solution of the homogeneous equation, u h9 and a 
particular solution of the non-homogeneous equation, w p , and hence 

" = «h+ w p (5) 

The particular solution « p can be written as 

00 

Wp 555 2 a n %mh n z (6) 

71=1 
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where a n is a complex constant. Substitution of equation (6) into equation (4) yields 

00 00 00 

£ P ^2 a n h*smh n z-maj 2 J a n sinh n z = - 2 a jm^n sinA »z 



rc-l 



n=l 



Thus, the constant a 7l can be obtained as 



0) 



(8) 



The solution of the homogeneous equation can be written as 

w h = ^sin \z+B cos Az+Csinh Az+Z> cosh Az (9) 

where ,4, 5, C and Z) are integration constants determined by the boundary conditions at the ends of the pile 
and 

A-^/(W/£ p 7) (10) 
Then by equation (5), the pile amplitude is in general 



00 Oiu U 

u(z) = A sin Az-fi?cos Az+Csinh Az+Z) cosh Az— Y ~ r ,, n n 0 sin/Lz 



(11) 



According to equation (3) and equation (20) in Reference 1 the displacement of the soil layer at the pile can 



be written as 



U(z)= 2 U n siah n z 



n=l 



(12) 



Since a perfect connection between the pile and the soil layer is assumed, the horizontal displacement of the 
soil layer, equation (12), and the displacement of the pile, equation (1 1), must be identical so that 

00 Q 00 

/* sin Az cos Az+CsinhAz + Z> cosh Az- £ = _. h . n n g sinA w z= £ l/ n sin/Lz (13) 

n=1 h v in^-miD' n =i 

Expanding sinAz, cosAz, sinhAz and coshAz into a Fourier sine series of argument h n z> equation (13) 
becomes 



— JYltO J 



sinA n z== S V n %inh n z 



71=1 



(14) 



where 



2 C H \ 

Fin = 7? Sin ^ z S * n A» Z dz 

« Jo 
2 f H 

^n^T? cos Azsin A n zdz 
Jo 

F -If* 



sinh Azsin A n zdz 



cosh Az sin /*„zdz 



(15) 



Equation (14) can be solved for U n which yields 



_ ^F ln +JF 2n +Cf 3w + i)F 4n 



(16) 



Substituting t/ n into equation (11), the amplitude of the pile motion is 

« h „04F ln + J BF 2n + CF 3n + DF 4n ) . 



w = ^4 sin Az+jB cos Az+Csinh Az+Z) cosh Az— 2 



sinA n z (17) 
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Here, only the integration constants A, B, C and D remain undetermined. It is convenient to introduce the 
following dimensionless parameters: 



'"S' ' 0= "' h = Hh n = \>n{2n-\) 



H 



where the wave velocity ratio 



(18) 



(19) 



and /> p ,p = mass densities of the pile and the soil respectively. Then, equation (17), rearranged with respect 
to the integration constants, is also 

u = ^sin Az - y 2 / lw sin h n z + ^cos Az - Y 2 / 2 n sin /r n z^ £ + ^sinh fe-rS/ 3n sin h n z j C 



+ ^cosh Af - y 2 / 4n sin A n z j Z) 



(20) 



where 




<*hn 



'In 



1 2n 



[ 4n 



(21) 



and the Fourier coefficients, given by equation (15), are after evaluation 

2 C H 

= ~tj I sin AzsinA n zdz 

■"JO 



^AcosAsinA n _ \ . t 
2 — ro to forA^£ n 

^2_ A2 n 

for A = /L 



r 2n 



2 fH 

= — I cos Az 
"Jo 



sin/*„zdz 



/ A sin fl ft sin A h n \ T , 

^ 0 for A = /J n 

A cosh A sin 



*3n = T? f sinhAzsin// n zdz = 2^^ 

2 CH / 
F 4n = — J cosh Az sin A n zdz = 21 



+A 2 



Asin/i n sinhA 
1M 



K \ 



(22) 



In equation (21), a hn — a hn /(7r/x) = dimensionless resistance factor given by equation (31) in Reference 1 
and analysed in detail in the same reference. It characterizes the soil resistance to the motion of the pile in 
the nth wave mode of the soil layer. The total soil resistance is a sum of the contributions from all the modes. 
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With the displacement of the pile described by equation (20), the amplitudes of the angle of rotation iff, the 
bending moment M and the shear force S are obtained from the standard relationships, 



*(*)- 



du 

dz 



d^w 
dz 2 



d*w 
dz 3 



(23) 



All these quantities are listed in a matrix form in Table I. The summations in Table I extend from one to 
infinity but sufficient convergence is obtained with about a hundred terms. 

For any boundary condition of interest, the equations given in Table I yield a number of linear algebraic 
equations sufficient to establish the integration constants A, B> C and D, With these integration constants 
established, the horizontal displacement follows from equation (20) and the rotation and internal forces are 
obtained from equations (23). 

The cases in which the pile is excited by a horizontal force and moment applied at the head are particularly 
illustrative and useful. The steady response to a unit force (or moment) acting with variable frequency repre- 
sents the so-called compliance. Conversely, such end forces (moments) which produce a single unit displace- 
ment (rotation) of the pile head define the dynamic stiffness of the soil-pile system at a certain frequency. The 
real part of this stiffness represents the elastic restoring force ot the soil-pile system while the imaginary part 
of the stiffness characterizes the damping of the system. With the stiffness known, the response of any 
structure supported by piles can be predicted. 4 Hence, the stiffness of the soil-pile system is of primary 
interest. 

DYNAMIC STIFFNESS OF THE SOIL-PILE SYSTEM 
The boundary conditions at the lower end of the pile are for the pinned tip 

w(0) = 0, M(0) - 0 (24) 

and for the clamped tip 

w(0) = 0, 0(0) = 0 (25) 

The stiffnesses at the level of the pile head are defined by a unit horizontal translation and by a unit rotation 
of the pile head. These displacements define stiffnesses k 8w k Mu , k M ^ and k s ^ where the first subscript 



(+) k 



+ u<H) = I 




(-)k 




Figure 2. Generation of stiffnesses of soil-pile system at level of pile head 

indicates the horizontal force 5 or moment M and the second subscript identifies the unit displacement of the 
pile head (Figure 2). Hence, the boundary conditions at the head are for k Su and k Mu 

W (f=l) = l, >A(z = l) = 0 (26) 
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and for k Mf and k Sf 



u(z= 1) = 0, xl>(z = 1) = 1 



(27) 



Under the assumptions adopted k S p = k Mu . Substituting the four boundary conditions chosen into the 
pertinent equations of Table I yields four linear equations for the determination of the integration constants. 
The first two of these equations are obtained by equations (24) or (25) and are for the pinned tip 




(28) 



and for the clamped tip 



0 



L A- rsUi 



1 0 



3n 



in 




0 

o J 



(29) 



The remaining two equations are obtained by equations (26) or (27) and are 

sin A — ySAn sm K cos A" - Y'Ef 2n sin h n 
A cos A - Y^Jinfxn cos h n — A sin A - Y 'EJi n f in cos h n 

sinhA-KS/sn 

sin h n cosh A — ^STin sin h n 
A cosh A - YYJi n f Zn cos h n A sinh A - Yh n 2 fi n cos h 




-o 



for fc Su and k 



Mu 



for k m and *^ 



(30) 
(31) 



For the clamped tip, equation (29) yields B = —D; equation (28) yields for the pinned tip B = D = 0. 

After obtaining the integration constants A, B 9 C and D for the frequency of interest, the stiffnesses are 
calculated from Table I with z = 1 and the following substitutions; 



ksu = S(H) and integration constants for u(H) = 1 
k Mu = —M(H) and integration constants for u(H) — 1 
k M ^ = M(H) and integration constants for $(H) = 1 
ksp =» —S(H) an d integration constants for </<//) = 1 



(32) 



Some numerical results showing the variations of the pile stiffness with the main dimensionless parameters 
are given later herein. The application of the boundary conditions to other cases, including harmonic excita- 
tion of the pile head, is analogous. 



STATIC STIFFNESS OF THE SOIL-PILE SYSTEM 

The static case can be obtained from the above relationships as a limit case for «>(A)->0 after the integration 
constants have been evaluated. However, no simple explicit formulas can be obtained for to = 0 from the 
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dynamic solution. If the static case is of particular concern it appears more suitable to develop a special static 
solution by modifying the above approach. 

In the static case, the displacement is given by equation (4) in which the term containing a> vanishes and 
thus 



,d 4 w 



(33) 



£p 1 d? " ? i ahn Un Sin kn Z 
The particular solution of equation (33) and the complete solution of its homogeneous part are respectively 

<*hn U n e 



Therefore, the complete solution is 



- sin h n z 
u h = Az*+Bz 2 +Cz+D 



Following the same process as that in the dynamic case, U n for the static case can be obtained as 

AF ln +BF 2n +CF 3n +DFt n 



where 



u n = - 



'In 



1 2n 



'3n 



sin h n z 




(34) 
(35) 

(36) 
(37) 

(38) 



Finally, introducing the dimensionless parameters, the horizontal static displacement can be written as 
u = (z»- Y £/ ln sin* n f^+^- YZf 2n smh n z^B+(z- Y js/ 8B sinfl n z)c+(l- Y js/ 4n 



xsin/J n £| 



(39) 



where 



and 







f F ln /H* \ 






F*JH* 


it j 




irl 









(40) 



(41) 



F ln = 2(^sin/i n ) 
F 2n = 2(^smh n -^j 



= 2^sin^„ 



The static horizontal displacement of the pile, its rotation, bending moment and shear are given in Table II. 
The integration constants A, B, C and D can be obtained in the same way as in the dynamic case described 
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above. For the calculation of the stiffnesses, equations (26) and (27) again specify the boundary conditions. 
The procedure is given in more detail in Appendix I. 



PARAMETRIC STUDY AND EXAMPLES 

In this section, some numerical results are presented for a few typical values of the dimensionless parameters 
characterizing the soil-pile interaction. 

The static stiffness 

The dimensionless parameters controlling the static stiffness of the soil-pile system depend only on the 
wave velocity ratio v* the slenderness ratio H/r 0 and the Poisson ratio v. For a broad range of these ratios, 
the static stiffnesses k 8u , k 8 ^ and k M ^ were calculated for both clamped and pinned tips and are shown, in a 
dimensionless form, in Figure 3. It can be seen that with increasing wave velocity ratio, that is with increasing 
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stiffness of soil, the dimensionless stiffnesses of the pile become less dependent upon the slenderness ratio 
and the tip conditions. The effect of the tip condition is strong for very short piles and weak soils. 

In Figure 4, the distribution of static horizontal displacement along the pile length is shown for a typical set 
of input parameters and the end conditions indicated. With a slender pile, only the upper part of the pile 
undergoes significant displacements, while the lower part remains effectively silent. Therefore, the soil 
properties used in the analysis should be characteristic of the uppermost layer of the soil in practical applica- 
tions. 

^ = 100 

To 

v = 0.0014 
^ = 0,6 
v = 0.5 

/////// 

Figure 4. Vertical distribution of static horizontal displacement of pile 
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Figure 5. Comparison of static displacements of pile head with Poulos' 8 solution 

To check the accuracy of the static solution, the results were compared with the finite difference solution by 
Poulos. 6 Static displacements of the pile head due to a unit horizontal load were computed for different soil 
modulae and plotted in Figure 5. Free as well as fixed heads were considered. The present solution is shown 
in Figure 5 in dashed lines and the solution by Poulos (his Figure 9) is shown in full lines. The agreement 
between the two solutions is very good. 



Dynamic stiffness and damping 

The dynamic stiffness of the soil-pile system is a complex number. Its real part represents the true (in 
phase) stiffness and its imaginary (out of phase) part represents the damping. The damping is generated by 
energy dissipation due to propagation of elastic waves to infinity and by energy absorption through its 
conversion into heat. The latter part of the damping is a result of internal friction in soil modelled here as 
frequency independent hysteresis associated with both volumetric and shear strains. These two components 
of internal damping can be characterized by internal damping ratios D Y and D a . The numerical examples 
presented in this paragraph were calculated assuming D V — D B = D. 
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The dynamic complex stiffness of the soil-pile system depends on the same dimensionless parameters as the 
static stiffness and, furthermore, on the internal damping ratio, the mass ratio p = p/p p and the dimensionless 
frequency. The dimensionless frequency is particularly illustrative if presented as a ratio b 1 = A/Aq, where A 
is given by equation (18) and A 0 = the lowest value of A at which the stiffness of a free-standing pile without 
soil becomes zero. The dynamic stiffness can also be shown in a dimensionless form as a ratio. 

h' = £(dynamic)/fc(static) 

The variations of the complex dynamic stiffnesses k r with some of the dimensionless parameters of the 
problem are shown in Figures 6-8. 




Figure 6. Variation of dimensionless stiffness constant k' Su with frequency, slenderness and tip conditions 

It can be seen that the real stiffnesses can diminish to zero and attain negative values when the piles are 
stiff (short) and the soil is weak. The smaller, sharp minima in the real stiffness appear at the first resonance 
with the soil layer. They dip to zero in the absence of material damping and, conversely, almost disappear 
with higher values of material damping than those used. Hence, it may be observed that the variation of real 
stiffness with frequency can be very dramatic for stiff piles, very weak soil and small material damping and 
remarkably modest for slender piles, stiff soil and large material damping or, in general, for low dimensionless 
frequencies. In some cases, even an increase in real stiffness with frequency can be observed. 

The influence of material damping on dynamic stiffness is further illustrated by Figure 9. Material damping 
has little effect on real stiffness below the first resonance with the soil layer, is very important in the first 
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several resonance regions where very sharp minima can occur if it is absent, and reduces the real stiffness at 
higher frequencies. This reduction of stiffness would be much more significant if frequency-dependent material 
damping was used. 

The imaginary part of stiffness representing the damping of the soil-pile system shows variations with 
frequency that often strongly differ from the linear law. Only in those frequency regions where these variations 
are close to linear ones can the damping be represented by a frequency-independent coefficient of equivalent 
viscous (Voigt) damping. Fortunately, such a possibility appears to be most acceptable for horizontal 




Figure 7. Variation of dimensionless constant k' s $ with frequency, slenderness and tip conditions 

stiffness in horizontal translation, k Sui which is most important. Very little damping is generated below the 
first resonance with the soil layer, because no progressive wave is generated in that frequency region in the 
absence of material damping and, consequently, no geometric damping occurs. Above the first resonance, 
the total damping is mostly due to geometric damping with little contribution from material damping. In the 
absence of material damping, sharp minima in the imaginary part of stiffness are associated with the sharp 
minima in the real stiffness (Figure 9). However, the markedness of this minima decreases with increasing 
order of the resonances and increasing stiffness of the pile. With stiff (short) piles, only the first resonance of 
the layer is usually marked (Figure 6). 

Effect of Poissorts ratio 

Figures 6-8 were calculated for one value of Poisson's ratio, v = 0-5. The variation of complex stiffnesses 
with Poisson's ratio is shown for the horizontal vibration and the theoretical limit cases of v — 0 and 0-5 in 
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Figure 10. The dashed curves represent the stiffness for one value of Poisson's ratio normalized by its static 
value and multiplied by the static stiffness for the other value of Poisson's ratio. Thus, the two bracketed 
curves actually represent the normalized stiffnesses k' computed for Poisson's ratios 0 and 0-5 if their values 
at a f 0 = 0 are viewed as being equal to unity. While the true dynamic stiffnesses, k, strongly depend on 
Poisson's ratio, the normalized stiffnesses, k\ are almost independent of Poisson's ratio. Representation of 
the dynamic stiffness is facilitated by means of a unique dimensionless stiffness valid for all Poisson's ratios 
and the static stiffness pertinent to any particular value of Poisson's ratio. 
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Figure 8. Variation of dimensionless constant k' M $ with frequency, slenderness and tip conditions 



0.4 0.6 0.8 _ _ 1.0 
FREQUENCY b ( = X/X 0 



1.2 



1.4 



A different situation exists with vertical vibration, as exemplified by Figure 11 showing true vertical 
stiffness, k, (a) and its values normalized by the static stiffness (b). The true vertical dynamic stiffness is almost 
independent of Poisson's ratio for higher frequencies, but depends on v for low frequencies and the static 
case. On the other hand, the normalized stiffness, k\ depends on Poisson's ratio in the vertical vibration. (The 
relationships shown in Figure 1 1 were computed from formulas given in Reference 2.) 



Dynamic response 

Figure 12 shows the response of the pile head to harmonic excitation for two tip conditions and two soil 
stiffnesses. Each response curve can feature two types of resonant peaks. The smaller peaks appearing at low 
frequencies are due to resonance with the natural frequency of the soil layer. The major peaks evident in the 
shorter pile are due to resonance with the soil-pile system. With very weak soils (v = 0-002) and stiff piles, 
these major peaks are well pronounced and occur close to the natural frequency of the pile alone (X/Xq = 1). 
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With slender piles, these major peaks can be completely suppressed by large radiation (geometric) damping 
generated. The smaller peaks due to resonance with the soil layer can become less marked when material 
damping is higher. The possible disappearance of the major peaks can render dynamic testing of piles without 
a cap impractical. 



Comparison with other dynamic solutions 

For comparison with other dynamic approaches, solutions by Tajimi 3 and Novak 4 were used. Comparison 
with the solution by Novak 4 indicates reasonable agreement for slender piles and higher frequencies (Figure 
5 in Reference 4). Comparison with Tajimi's solution 3 is excellent for the vibration mode analysed by him,* 
if the difference in the definition of material damping is accounted for. 




Comparison of the theory with experiments would require a special study. Experiments conducted with 
model pile foundations in the field and described in Reference 7 indicate that the linear elastic theory can 
give very good results for small amplitudes when a shear modulus of soil corresponding to the upper layer of 
soil is used, and when a correction for the effect of pile grouping (interaction) is made. An apparent modulus 
of soil derived from a static test was found useful as it can approximately account for the effects of pile 
separation from the soil and for the diminishing of the soil modulus towards the ground surface due to 
reduction of the confining pressure. 



CONCLUSIONS 

Interaction between soil and an elastic pile in horizontal vibration was theoretically examined. The solution 
was brought to a form which makes it possible to predict the dynamic displacements, internal forces, stiffnesses 
and damping for different soils and end conditions of the pile. A special solution for the static case was also 



Tajimi assumes horizontal motion of a pile with a clamped head and a pinned tip and viscous damping. 
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developed. The computing costs of numerical calculations are quite modest.* The main observations made can 
be summarized as follows: 

A direct solution of the response appears superior over the modal analysis of forced oscillations. 

The solution presented agrees well with other solutions available for both the static and dynamic cases. 

The analytical solution offers an interesting insight into the complicated nature of soil-pile interaction. 




Figure 10. Effect of Poisson's ratio on complex stiffness of pile in horizontal vibrations 

The dynamic stiffness and damping of the soil-pile system depends on a number of dimensionless parameters 
such as the slenderness ratio, the wave velocity ratio, dimensionless frequency, Poisson's ratio, the internal 
damping ratio and, to a lesser degree, the mass ratio. 



* The computer program can be obtained from SACDA, the Systems Analysis, Control and Design Activity, Faculty of 
Engineering Science, University of Western Ontario, London, Ontario, Canada. 
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Under some conditions both the stiffness and damping exhibit dramatic variations with frequency. The 
dynamic stiffness can diminish to zero and even achieve negative values. 

The effect of material damping is very significant in the region of resonance with the lowest natural 
frequencies of the soil layer. Above this region, small material damping has little significance but large 




Figure 12. Response of pile head to horizontal excitation and corresponding phase lag 

material damping reduces the stiffness progressively with increasing frequency; however, the total damping is 
produced primarily by geometric damping. Below the first natural frequency of the layer, all damping gener- 
ated is due to material damping and hence, no damping occurs in absence of material damping. 
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APPENDIX I 

Calculation of the integration constants for static stiffness 

The first two equations for the determination of the integration constants are for the pinned tip 

0 0 0 1 1 5 (0 
0 2 0 0 \ C I 0 



(42) 



and for the clamped tip 



0 



/In 



0 0 
- Y^hn/zn 1 - Y^fl n f 3n 



L - Y^hJ, 

The remaining two equations are by equations (26) or (27) 




(43) 
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for k Su and k- 



Mu 



(44) 



for k Mf and k Sf 



With the integration constants established the static stiffness follows from Table II by equations (32). (It 
may be mentioned that the integration constants are not dimensionless.) 



APPENDIX II 



Notation 



A,B,C,D integration constants 

a n complex constant of the soil layer in horizontal direction 

aQ — r 0 co/v B dimensionless frequency 

b 1 = X/Xq ratio of excitation frequency to frequency Ao 

D assumed damping ratio for D Y = D a 

D y , D 8 damping ratios associated with volumetric and shear strain respectively (Z) v = A'/ A, 

A = /*'//*) 

E Young's modulus of soil 

E v Young's modulus of pile 

f in dimensionless coefficient 

F in Fourier coefficient 

H length of pile, depth of soil layer over the bedrock 

h n (7772//) (2«-l), n = 1,2,... parameter 

h n — Hh n dimensionless parameter 

/ second moment of area of pile cross-section 

i = ^/ — 1 imaginary unit 

ky stiffness in direction i corresponding to unit displacement in direction j 

k! = &(dyn)//;(stat) dimensionless dynamic stiffness 
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m mass of pile per unit length 

p(z) amplitude of resistance force of the soil layer 

r 0 radius of pile 

r 0 = rJH ratio of pile radius to layer thickness 

5 shear force, (cross-sectional area of pile in Figure 11) 

t time 

U n modal amplitude of pile displacement 

u = u(z) complex amplitude of pile displacement 

£ = v sl v p wave velocity ratio 

v p = V(£p//° P ) longitudinal wave velocity in pile 

v s = <J(jjl/p) shear wave velocity in soil 

Y = 4v* pjrft dimensionless parameter 

z vertical co-ordinate 

z = z/H relative co-ordinate (depth) 

<x hn resistance factor of a soil layer in horizontal direction 

&hn ~ 0L hJ 7T p< dimensionless resistance factor in horizontal direction 

A, /x real part of Lame's constants for soil 

A', p! imaginary part of Lame's constants for soil 

^ = ^[(/kw 2 /£ p /)] frequency parameter 

A = HX dimensionless frequency 

Ao dimensionless frequency at which stiffness of a free-standing pile becomes zero 

fx shear modulus of soil 

v Poisson's ratio of soil 

p p mass density of pile 

p mass density of soil 

p = p/p v mass ratio 

a> circular frequency 
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